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Abstract 

00 

The two-loop self-energy correction is evaluated to all orders in Za for the ground- 
state Lamb shift of H-like ions with Z > 20, where Z is the nuclear charge number 
and a is the fine structure constant. The results obtained are compared with ana- 
£^ ■ lytical calculations performed within the Za-expansion. 
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^ ■ The subject of this paper is a study of the historically most problematic two- 

loop QED effect for the Lamb shift of H-like ions, the two-loop self-energy 
correction. Traditionally, investigations of radiative corrections to a given or- 
der in the fine-structure constant a rely on an expansion in the parameters 
Za and ln[(Za)~ 2 ] (Z is the nuclear charge number). For the two-loop self- 
energy correction, the energy shift is conveniently expressed in terms of a 
dimensionless function F(Za), 



AE = m (*) 2{ -^F(Za), 
\7r/ n 6 

which expansion reads 

F(Za) = B m + (Za)B 50 + {Za) 2 \L 3 B 63 + L 2 B m + L B 61 + G Ko (Za) 



Email address: yerokhin@pcqntl.phys.spbu.ru (V. A. Yerokhin). 



Preprint submitted to Elsevier Science 



2 February 2008 



(2) 



where L = ln[(Za)~ 2 ], n is the principal quantum number, and the remainder 
G h , XZa) incorporates all higher-order terms. Calculations of the coefficients 
of this expansion extended over several decades [1,2,3,4,5,6], 



B 40 (ns) = 1.409 244... , (3) 

B 50 (ns) = -24.2668(31), (4) 

£? 63 M = -8/27, (5) 

fl 62 (ls) = 16/27 - (16/9) In 2, (6) 

B 61 (ls) =49.838317. (7) 

The expansion of the higher-order remainder Gh.o. starts with a constant, 
G^. .{Za) = BQ +Za(- • •) . A partial result for the B 6Q coefficient was obtained 
recently by Pachucki and Jentschura [7], 

B m = -61.6(3) ±15%, (8) 



where ±15% refer to the possible contribution of uncalculated terms. One of 
the notable features of the expansion (2) is its remarkably slow convergence, 
even for low values of the nuclear charge number Z, which exhibits itself by 
large numerical values of the expansion coefficients B^. This makes clear that 
a direct evaluation of the two-loop self-energy correction is highly desirable, 
especially when addressing middle- and high-Z ions. 

Calculations of the two-loop self-energy correction (Fig. 1) to all orders in 
the parameter Za started with the irreducible contribution of the diagram 
(a) (also known as the loop- after-loop (LAL) correction), which is by far 
the simplest part of the total set. Such evaluation was first accomplished 
by Mitrushenkov et al. [8] and later by other authors [9,10]. The contribu- 
tion of the remaining diagrams in Fig. 1 is by far more difficult to calculate. 
First attempt to evaluate it to all orders in Za was made by Mallampalli and 
Sapirstein [11]. In that work, the contribution of interest was rearranged in 3 
parts, referred to by the authors as the "M", "P", and "F" terms. Mallam- 
palli and Sapirstein calculated only the M and F terms, while the P term 
was left out since a new numerical technique had to be developed for its com- 
putation. In addition, since the numerical procedure turned out to be very 
time consuming, the actual calculation of the M term was carried out only 
for two ions, uranium and bismuth. Subsequently, in the investigation by two 
of us [12], a computation of the remaining P term for Z = 83, 90, and 92 
was performed, which formally accomplished the calculation of the two-loop 
self-energy for these ions. 

First complete evaluation of all contributions to the two-loop self-energy cor- 
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rection to all orders in Za was carried out in our recent investigations [13,14] 
for the ground-state Lamb shift of H-like ions with Z > 40. In the present 
work, we extend our calculation to the ions with Z = 20, 30 and improve the 
accuracy of the previous numerical results for Z = 40-60. 

It is beyond the scope of the present paper to describe all modifications of the 
numerical procedure that allowed us to obtain improved numerical results in 
the low-Z region, despite the occurrence of large numerical cancellations that 
tend to grow very fast as Z decreases. We only mention that the new dual- 
kinetic-balance basis set [15] constructed with B splines was employed in our 
calculations of the P term, which considerably improved the convergence with 
respect to the number of basis functions in the low-Z region. Another problem 
that was encountered in the evaluation of the P term was to keep numerical 
integrations over momentum variables well under control. This problem is 
associated with a significant contribution coming from the region of very large 
momenta, where the numerical Green function is not smooth enough, due 
to restrictions of a finite-basis-set representation. The problem was solved 
by introducing a set of subtractions that have the same behaviour for large 
momenta as the original integrand but are easier to evaluate numerically, 
and by employing very fine grids for numerical momentum integrations. The 
calculation of the M term was carried out employing the contour Clh for the 
integrations over the virtual photon energies (see Ref. [14]) that is much more 
suitable for the numerical evaluation in the low-Z region than the integration 
simply along the imaginary axis. In addition, it turned out to be possible to 
separate from the M term a contribution, which contains the dominant part of 
the spurious terms in the \ow-Z region that are cancelled in the final sum, and 
to calculate it separately. This contribution contains only one free partial- wave 
expansion parameter, which makes the corresponding numerical computation 
much easier as compared to the full M term that has a partial-wave expansion 
over two independent parameters. 

Numerical results obtained for the two-loop self-energy correction for the 
ground state of H-like ions are presented in Table 1 in terms of the func- 
tion FiZa) defined by Eq. (1). As compared to our previous investigations 
[13,14], we performed calculations for Z = 20, 30 and improved the numerical 
accuracy of the results for the P and M terms for Z = 40, 50, and 60. The 
values for Z > 70 are identical to those in our previous evaluation [14] . 

In Fig. 2 we compare our non-perturbative (in Za) results with contributions 
of the known coefficients of the expansion (2). The solid line corresponds to the 
first two terms in Eq. (2), whereas the dashed line represents a contribution of 
all known terms of the Za expansion (including the result of Ref. [7] for Bqo). 
In order to provide a more detailed comparison with the analytical results, 
we first separate the contribution containing all orders in Za starting from 
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(Za) 2 , defined as 



F(Za) = r ^- n *° =B 50 + (Za)(...). (9) 



The function FiZa) is plotted in Fig. 3. The cross on the picture represents 
its analytical value at Z — 0, B 50 . It is noteworth that, despite the numerically 
large contributions Bg 3 -B m in the next-to-leading order, the behaviour of the 
function F(Za) is very smooth and that it tends to the analytical value at 
Z = 0. 

Next, we separate the higher-order remainder G h '°'(Za) [defined by Eq. (2)] 
from our non-perturbative results, with the corresponding plot presented in 
Fig. 4. The cross with the error bars corresponds to the analytical result of 
Ref. [7] for the Bqq coefficient, Bqq = —62 ±9. We see that the apparent limit 
of the numerical values at Z — >■ seems to be about twice as large as the 
analytical result of Ref. [7]. However, we can not at present conclude whether 
this disagreement stems from the coefficient B 60 or from one of the logarith- 
mic coefficients in the order a 2 (Za) 6 . A detailed comparison with individual 
coefficients can be in principle performed by fitting numerical results to the 
analytical form of the Za expansion, but it requires an improvement of the ac- 
curacy of the numerical data. We mention also that a reliable extrapolation of 
the numerical data to Z = is complicated by the presence of nearly-singular 
terms like ZalnZa, Zaln 2 Za in the Za expansion of the remainder G ho ', 
which could possibly result in a rapidly-varying structure of G h o (Z) near the 
origin. We hope that an extension of our calculations to lower values of Z and 
an improvement of numerical accuracy of the results will allow us to provide 
a more detailed comparison with analytical calculations. 

In summary, we have performed calculations of the two-loop self-energy correc- 
tion that extend our previous evaluations to the region Z > 20. The numerical 
results obtained are compared with the known terms of the Za expansion. A 
good agreement is demonstrated with the first terms of this expansion. How- 
ever, the numerical results are shown to be likely in a disagreement with the 
analytical results to order a 2 (Za) 6 . 

Valuable discussions with K. Pachucki and U. Jentschura are gratefully ac- 
knowledged. This work was supported by INTAS YS grant No. 03-55-1442, 
by RFBR grant No. 04-02-17574, and by foundation "Dynasty". The com- 
putation was partly performed on the CINES and IDRIS national computer 
centers. 
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(a) (b) (c) 

Fig. 1. Two- loop self-energy diagrams. 
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Table 1 

Individual contributions to the two-loop self-energy correction expressed in terms 
of F(Za). 



z 


LAL 


F term 


P term 


M term 


Total 


20 


-0.602 


136.89(3) 


-102.03(9) 


-34.77(9) 


-0.52(13) 


30 


-0.757 


44.723(5) 


-29.42(4) 


-15.46(4) 


-0.91(6) 


40 


-0.871 


19.504(2) 


-11.58(2) 


-8.26(7) 


-1.21(7) 








-11.41(15) a 


-8.27(18) a 


-1.05(23) a 


50 


-0.973 


10.025(2) 


-5.486(20) 


-5.00(5) 


-1.43(5) 








-5.41(8) a 


-4.99(6) a 


-1.34(10) a 


60 


-1.082 


5.723(1) 


-2.961(20) 


-3.342(9) 


-1.66(2) 








-2.93(4) a 


-3.342(21) a 


-1.63(4) a 


70 


-1.216 


3.497 


-1.757(25) 


-2.412(11) 


-1.888(27) 


83 


-1.466 


1.938 


-1.057(13) 


-1.764(4) 


-2.349(14) 


92 


-1.734 


1.276 


-0.812(10) 


-1.513(3) 


-2.783(10) 


100 


-2.099 


0.825 


-0.723(7) 


-1.384(3) 


-3.381(8) 



Ref. [14]. 
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Fig. 2. The two-loop self-energy correction to all orders in Za (dots) together with 
the contribution of the first two terms of the Za expansion (2) (solid line) and all 
known terms of the Za expansion (dashed line). 
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Fig. 3. Numerical results for the two-loop self-energy correction with the leading 
contribution of the Za expansion subtracted out, as given by Eq. (9). The cross 
denotes the analytical value for this contribution at Z = (i.e., the coefficient i?5o)- 




Fig. 4. Higher-order remainder G La (Za) defined by Eq. (2). The cross denotes the 
analytical result for this contribution at Z = (i.e., the coefficient Bqq) obtained in 
Ref. [7]. 
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